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Abstract 

Applying the Amplitude Difference method to Fermilab experiment E791 
D + — ► ir~Tc + ir + data, we measure the low mass tt + tt~ phase motion. Our 
results suggest a significant phase variation, compatible with the existence of 
an isoscalar <r(500) meson, as previously reported using an isobar model fit to 
the full Dalitz-plot density. 
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Recently we proposed the isobar-based Amplitude Difference (AD) method to 
extract the phase motion of a complex amplitude in three-body heavy-meson decays 
[1]. With this method, the phase variation of a generic complex amplitude can be 
directly revealed through interference in the Dalitz-plot region where it crosses a well 
established resonant state, used as a probe. As a test, this method was successfully 
applied to data [2] to extract the well known phase motion of the scalar amplitude 
/o(980) observed in D+ — > it~it + it + 1 decay. In the present paper we use the same 
method to study the low 7r + 7r~ mass region of the D + — > 7r~7r + 7r + decay where 
Fermilab experiment E791 showed evidence for the existence of a light and broad 
scalar resonance [3]. 

To obtain good fit quality in a full Dalitz-plot analysis, E791 found it necessary to 
include an extra scalar particle, in addition to the well-established di-pion resonances 
[4]. For this new scalar state, parameterized as an S-wave Breit-Wigner resonance, 
they measured a mass and a width of 478^23 ± 17 MeV/c 2 and 324±^ ± 21 MeV/c 2 , 
respectively. These parameters are compatible with those expected for the isoscalar 
meson <j(500). The D + — > <r(500)7T + [3] decay appeared as the dominant contribution, 
accounting for approximately half the D + — > 7r~7r + 7r + decays. 

The E791 result has been widely discussed [5]-[10] and new data has become 
available [11]. However it is desirable to be able to confirm the result through a 
direct observation of the phase motion expected for a resonance [6, 7, 8]. In this 
context we apply the AD method to the low 7r + 7r~ mass region of the D + — > 7r~7r + 7r + 
decay. We also compare the phase variation of the Breit-Wigner function found in 
the isobar Dalitz-plot analysis [3] to the model-independent method of this paper. 

The present study is a reanalysis of the Fermilab experiment E791 data. Here 
we investigate a subset of the total phase space used by the experiment in their 
full Dalitz-plot analysis. A description of the experiment, data selection criteria, 
background parametrization and detector acceptance are found in references [3, 12]. 
The final 7r~7r + 7r + invariant mass distribution is shown in Fig. 1. There are 1686 
events with invariant mass between 1.85 and 1.89 GeV/c 2 shown in the shaded region 
of Fig. 1. The integrated signal-to-background ratio in this range is about 2:1. Fig. 2 
shows the folded Dalitz-plot. The horizontal and vertical axes are the squares of the 
7r + 7r~ invariant mass high (s 12 ) and low (s 13 ) combinations. The analysis presented 
here uses the hatched area of Fig. 2. We estimate 60 background events in a total 
of 197 candidate events. The background does not show any dependence on the Si 2 
variable [3]. 

The detector acceptance in this region is almost constant. There is a very mild 
slope in the S13 acceptance, but no variation with si 2 . Nevertheless the acceptance is 
taken into account to correct the event distributions shown later. 

There are two conditions necessary to extract the phase motion of a generic am- 

1 Charge conjugate states are implied throughout the paper. 
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Figure 1: The n~n + n + invariant mass spectrum. The dashed line represents the total 
background. Events used for the D + isobar Dalitz-plot analysis [3] and in this AD 
method analysis are in the hatched area. 

plitude with the AD method: 

• A crossing region between the amplitude under study and a probe resonance 
has to be dominated by these two contributions. 

• The integrated amplitude of the probe resonance must be symmetric with re- 
spect to an effective mass squared (m^). 2 

To study the low mass region in Si 3 , three well known resonances could serve as 
a probe in s 12 in the D + -> decay: p(770), /o(980) and / 2 (1270). How- 

ever the broad p(770) and /o(980) are too close to each other to pass the isolation 
criteria mentioned above, as can be seen in Fig. 2. On the other hand, the tensor 
/ 2 (1270), m\ = 1.61 GeV 2 /c 4 , is located where the p(770) reaches a minimum due 
to its decay angular distribution in the crossed (si 3 ) channel, (see Fig. 3). In the 
D + — > 7t~7t + 7t + decay, the /2(1270) contribution satisfies the necessary conditions of 
having a substantial contribution crossing the low mass region, this being the regions 
where all other amplitudes can be considered negligible. In particular, we estimate a 
contamination of 5% of p(770)yr + events in the region of interest. 

We assume that the only contributions in this region are the / 2 (1270) amplitude 
in S12 and the mi complex amplitude under study in si 3 . We write: 

A(S 12 ,S 13 ) « a R £W/ 2 (1270)(S12) J=2 -A/f/ 2 (l270)(si2,Sl 3 ) + (1) 

+_a s /{p*/^- Sl ,) sinS(s 13 ) 

2 An analysis with more amplitudes in a limited phase space could be done, but in a model 
dependent way. 
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Figure 2: The folded D + — > 7r~7r + 7r + Dalitz-plot distribution of the events in the 
hatched area of Figure 1. The events used by the AD method analysis are in the 
hatched region of this folded Dalitz-plot. The size of the area of each bin in the plot 
corresponds to the number of events in that bin. 



where 7 is the overall relative phase-difference between the two amplitudes (assumed 
to be constant and arising from both production and final-state-interaction (FSI) be- 
tween the di-pion system and the bachelor pion); sinS(si^)e lS ^ Sl3 ^ represents the most 
general amplitude for a two-body elastic scattering, p*/\^s 13 is a phase space integra- 
tion factor to make this description compatible with mr scattering and p* is the pion 
momentum measured in the resonance rest frame; J=2 .M/ 2 (i27o)(si2, S13) is the angu- 
lar function for the / 2 (1270) tensor resonance given by |(| p 3 || p 2 \) 2 (3cos 2 9 — 1), 
9 being the angle between the pions 2 and 3, and J the angular momentum of the 
resonance; and clr and a s are the production strengths. The AD method assumes 
that the 7r + 7r~ production is constant over the effective mass range considered and 
clr and a s are energy independent. Finally the Breit-Wigner distribution is given by: 



BW = , ( 2 ) 



m T 

wiq — s — im r(s) 

,2,7+1 j 1 



The width is given by T(s) = T ^ jf 2 (^*) ' wnere the central barrier 

factor is J=2 F = 1/^9 + 3(rp*) 2 + (rp*) A . The parameter r is the radius of the 
resonance (~ 3 /to) [13] and p* = p*{m) is the momentum of decay particles at mass 
to, measured in the resonance rest frame, Pq = p*(m ), where to is the resonance 



mass. 
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Both r(s) and the angular function J=2 Mf 2 (i27o) for the / 2 (1270) resonance pro- 
duce asymmetries in the s l2 distribution. We use a D + — > / 2 (1270)7r + fast Monte 
Carlo (MC) simulation 3 to study the behavior of the probe resonance distribution 
in the Si 2 and S13 variables. The event distribution for / 2 (1270) in the s 12 subsys- 
tem from Monte Carlo is shown in Fig. 4. We choose an effective square mass of 
1.535 GeV 2 /c 4 such that the number of events between m 2 ^ and m 2 ^ + e (e = 0.26 
GeV 2 /c 4 ) is equal to the number of events integrated between m 2 ^ and m 2 ^ — e. 




Figure 4: Fast MC simulation/^ (1270) s 12 distribution, divided into m^ + t and m^-e, 
integrated between the threshold and 0.8 GeV 2 in s 13 . 

3 The Monte Carlo events are generated according to the physical amplitude squared using a 
uniform phase space density and weighted by the detector acceptance distribution over the Dalitz- 
plot. 
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For events in si 2 between m^j and fn 2 ^ + e in si 2 , J=2 A / l/ 2 (i27o)(si3) is shown 
in Fig. 5a and for those between m 2 ^ and rn 2 ^ — e in Fig. 5b 4 . We can see 
that these two plots are just slightly different. In our analysis we consider the ap- 
proximation J=2 -M- J f 2 (\27o) ( s i3) ~ J=2 -^/ 2 (i270)( Sl 3) an< ^ take ^he avera g e function 
J=2 7W/ 2 (i27o)(si3)- An important effect that we have to take into account is the zero 
of this function at s\ 3 ~0.48 GeV 2 /c 4 . Below we discuss the consequences of that for 
this AD method application. 




Figure 5: Fast MC of angular distribution function alone, J 2 'M. / 2 (i27o) , i n s i3 f° r ( a ) 
events between m 2 ^ and rn 2 ^ + e and (b) events between m 2 fj and m 2 ^ — e. 

To be brief, from here on we use J=2 -M/ 2 (i27o)(si 3 ) = M- and p*/y/s 13 = p'. 
The main equation of the AD method applied to the integrated amplitude-square 
difference is [1, 2]: 

/_ rm 2 . f +e fm 2 ,, 
A 2 ) = p'/M( / I A(s 12 , s 13 ) | 2 ds 12 - / ; n I A(s 12 , s 13 ) | 2 ds 12 ) 

Jm eff Jm efl- e 

~ -C(sin(25(s 13 ) + 7) - siwy), (3 

where C is an overall constant. 

From Eq. 2, it follows that the function A / A 2 p' / M. directly reflects the behavior 
of 5(si S ). A constant A / A 2 p' / M implies a constant S(s 13 ) which is the case for 
a non-resonant contribution. In the same way, a slow phase motion will produce a 
slowly varying A / A 2 p'/ M. distribution, and a full resonance phase motion produces 
a clear signature in A / A 2 p' / Ai with the presence of zero, maximum and minimum 
values. 

The plots were produced by a fast MC of the angular distribution function alone, with no 
detector influence included. We represent this function with the same binning used for our data 
events. 
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As mentioned previously, the background has no dependence on s\ 2 and its contri- 
bution vanishes from the A J A 2 distribution. The / A 2 in s 13 for events integrated in 
S12 miff and m 2 e ^ + e and m 2 e ff and miff — e, corrected by the acceptance shape, are 
presented in Figs. 6a and 6b, respectively; these events correspond to the hatched 
area of Fig. 2. We obtain A J A 2 shown in histogram 6c, by subtracting the 6b 
histogram from that in 6a. 




0.25 



Figure 6: Event distributions projected onto the S13 axis (a) for all events in the 
S12 interval J™£ +e | A(s 12 ,s 13 ) \ 2 ds 12 , and (b) for events f™ 2 °_ e | A(s 12 ,s 13 ) \ 2 ds 12 . 
The distributions are acceptance corrected such that the overall data statistics is 
conserved. Plot (c) shows the A / A 2 distribution . 

To extract the phase motion, we divide A J A 2 (Fig. 6c) by Ai (average of the 
distributions in Figs. 5a and b), and multiply by p', both known functions of s\ 3 . 
Then the only s± 3 dependence of the right hand side of Eq. 2 is in the phase motion 
S(sis). From Fig. 5, the zero at s 13 ~ 0.48 GeV 2 /c 4 in the angular function, produces 
a singularity around this value in A / A 2 p'/M. In Fig. 7 we show the A / A 2 p' / 
M. distribution. To treat the effect of the singularity, we have used a binning such 
that the singularity is placed in the middle of one bin. Doing this, we isolate the 
singularity in a single bin (bin 6) and discard it in the analysis. We point out that 
the location of this singularity can only affect the exact position of the minimum of 
A / A 2 p'/M. It does not change the general features observed, in Fig. 7, that this 
quantity starts at zero, has maximum and minimum values, and comes back to zero, 
the signature for a strong phase variation. The confidence level for a straight line 
fit to the data in Fig. 7 is 4.6%, while the separation between the maximum and 
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minimum values has a significance level of 2.6 r.m.s. 

We can see that the 6 th bin has a huge error, which corresponds to the bin due to 
the presence of the singularity. 
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Figure 7: The distribution of A / A 2 p'/A4 in s 13 . Note that the scale is arbitrary 
once the data, A j A 2 , is multiplied by an unnormalized function p' / Ai. The error 
bars represent statistical errors only. 

Assuming that £(513) is an analytical function of S13, Eq. 2 allows us to set 
the two following conditions at the maximum and minimum values of A / A 2 p'/A4, 
respectively. 

(A J A 2 )p'/M max -> sin(25(s 13 ) + 7) = -1 (4) 

and 

(A J A 2 )p'/M mm - S m(2<5( Sl3 ) + 7) = 1. (5) 

With these two conditions, we obtain 7 and C, and with these values we can extract 
directly 5(si3) from Fig. 7 by inverting Eq. 2. 

To propagate the statistical errors from Fig. 6 to the values of the 7 and 5(sis), we 
"produce" statistically compatible "experiments" by allowing each bin of Fig. 6a and 
Fig. 6b to fluctuate randomly following a Poisson law. We then solve the problem 
for each set. The statistical error in each bin for 5(si 3 ) is the r.m.s. of the S(s 13 ) dis- 
tributions from the Monte Carlo experiments. For the systematic errors, we change 
the e parameter (e = 0.22 GeV 2 /c 4 and e = 0.30 GeV 2 /c 4 ); we examine the possi- 
ble influence of other neglected amplitudes contributing in this region of the phase 
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space (based on the E791 amplitude measurements for non- resonant, /o(1370)7r + and 
p°(1450)7T + [3] contributions); and, to study the effect of averaging J=2 .M/ 2 (i27o)(si3) 
distributions, we use each distribution of Fig. 5a and Fig. 5b separately. The three 
systematic errors, while treated separately bin-by-bin, are found to be of an average 
size, relative to the statistical uncertainty, of 1, 0.6, and 0.4, respectively. They are 
added in quadrature. 

We measure 7 = 3.31 ± 0.33 ± 0.49 (with the first error statistical and the second 
systematic) 5 . The value is somewhat larger than the E791 full Dalitz-plot analysis 
value {^Daiitz = 2.59 ±0.19) [3]. The asymmetry of the distribution in Fig. 4 and the 
consequent use of an effective mass-squared for the ^2(1270) = 1.535 GeV 2 /c 4 instead 
of the nominal mass is responsible for the observed shift. To evaluate the magnitude 
of this effect, we generated 1000 fast MC samples with only two amplitudes, / 2 (1270) 
and <r(500). For both, we used Breit-Wigner functions with the E791 parameters, 
including the phase- difference of 2.59 rad. We extract 7 from these 1000 samples 
with the method presented here. The result has a mean value of 3.07 ± 0.10 rad, 
instead of the input value 2.59 rad. We estimate an offset of -0.48 (2.59 - 3.07) for 7 
from the difference between the generated and measured values in this Monte Carlo 
test. This yields to a corrected 7 cor r — 2.83 ± 0.38 ± 0.49. The production phase- 
difference between the /2(1270)7r + and the <t(500)7t + decays of D + measured in the 
isobar Dalitz analysis is in good agreement with 7 cor r from the AD method. 

With our 7 and C values we solve Eq. 2 for 5(s 13 ) for each si 3 bin. However, 
there are ambiguities that arise due to the sin" 1 operations. Table 1 shows the four 
possible solutions for 5(si 3 ) 6 . To resolve the ambiguities we use the assumption that 
the phase is zero at threshold and is an increasing, monotonic, and smooth function 

Of S13. 

The solution, after the above criteria for (5(si3), including systematic and statistical 
errors, is shown in Fig. 8 and in bold values in Table 1. We see a strong phase variation 
of about 180°, starting at threshold and saturating around S13 = 0.6 GeV 2 /c 4 . The 
limited sample size does not allow us to perform an accurate measurement of the 
mass and width parameters of a BW resonance fitted to this plot. We can only 
say that, at least for the preferred solution, the phase variation is compatible with 
the complete phase motion through 180° expected for a resonance. This result is in 

5 Here, in addition to the systematic error sources mentioned above, we have also incorporated 
one related to the binning choice that we made to avoid the singularity in A4. If we change the 
binning such that the singularity is placed between bins we measure 7 = 3.52 ± 0.53. 

6 The occurrence of two different sets of statistical errors in <5(si3) is caused by convolution of the 
error in 7 and in A J A 2 , and the fact that the distributions of those variables are not symmetrical. 
Moreover, the method naturally favors <5(si3) values coming from the maximum and minimum bins 
in Figure 7, since these are used to determine the C and 7 that are in turn used to determine all 
<5(si3) values. The uncertainties are therefore relatively smaller in these two bins than in the other 
bins, a result that may be attributed in part to the model itself rather than the quality of the data 
in these bins. 
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Table 1: Second column are the A J A 2 p' / M. values plotted in Fig. 7; note the arbi- 
trary normalization. Columns 3 to 6 are four possible solutions for 5(si3): <5o( s i3) for 
(25(s 13 ) + 7), <5i(si 3 ) for (2<J(si 3 ) + 7) + 27r > ^(sis) for vr - (25(s 13 ) + 7)) and 5 3 (s 13 ) 
for 7r — (25(s 13 ) + 7)) + 2ir). The systematic errors are described in the text. 



agreement with the evidence for a broad and low mass scalar resonance suggested by 
the previous E791 result using the full Dalitz-plot analysis [3]. 

The Breit-Wigner phase motion for the E791 mass and width parameters, com- 
bined with other scalar contributions obtained from that fit [3], is shown as the 
continuous line of Fig. 8. We can see a two standard deviation difference at lowest 
and highest 7r + 7r~ invariant masses squared, in opposite directions, but in overall 
agreement with the mass region where a scalar meson a must have its strong phase 
variation. Both results, the E791 with a Breit-Wigner phase variation, and the one 
presented in this paper, show a stronger phase variation than that obtained with the- 
oretical constraints in irn — > nn elastic scattering data in the scalar-isoscalar channel 
below 1 GeV [9, 14, 15]. The discrepancy between these results could be an indica- 
tion that applying Watson's theorem [16] is not straight-forward when comparing the 
phase motion of a two-body elastic interaction to the three-body decays. 

We have presented an extraction of the phase motion of the low mass tt + tt~ scalar 
amplitude using the well known /2(1270) tensor meson in the crossing channel acting 
as an interferometer. The result is obtained with an event counting procedure in a 
region of the phase space which is dominated by a D-wave / 2 (1270) interfering with the 
S-wave. The derivation of the phase motion relies heavily on the assumption that the 
maximum and minimum bins in Fig. 7 correspond to the quantity S = sin(25 ± 7) = 
+ 1 and —1. The clear presence of a maximum and a minimum, separated from each 
other by 2.6 standard deviations, supports this assumption. Given this caveat, the 
solution for <5(si 3 ) has a variation of about 180°, consistent with a resonant (x(500) 
contribution. We also obtain good agreement between the FSI 7 corr determined with 
the AD method and the 7 observed in the full Dalitz-plot analysis using an isobar 
model [3]. These results support the previous evidence for an important contribution 
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Figure 8: The phase values <5(si3) from our preferred solution versus the invariant 
7r + 7r~ mass squared with statistical and systematic errors added in quadrature. The 
continuous line is the Breit-Wigner phase motion with the E791 parameters for the 
(7(500) [3]. 

of the isoscalar <r(500) meson in the D + — > 7r~7r + 7r + decay [3]. 
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